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Abstract 

The inactivation of a Na+ channel occurs when the activation of 
the charged S4 segment of domain DIV is followed by the binding of 
an intracellular hydrophobic motif which blocks conduction through 
the ion pore. The voltage dependence of Na+ channel inactivation 
is, in general, dependent on the rate functions of the S4 sensors of 
each of the domains DI to DIV. If the activation of a single volt¬ 
age sensor that regulates the Na-I- channel conductance is coupled 
to a two-stage inactivation process, the rate functions for inactiva¬ 
tion and recovery from inactivation, as well as the time dependence 
of the Na-|- current in terms of the variables m(t) and h(t), may 
be derived from a solution to the master equation for interdepen¬ 
dent activation and inactivation. The rate functions have a voltage 
dependence that is consistent with the Hodgkin-Huxley empirically 
determined expressions, and exhibit saturation for both depolarized 
and hyperpolarized clamp potentials. 
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INTRODUCTION 

The opening and subsequent inactivation of Na+ channels and the activation 
of K+ channels generate the action potential in nerve and muscle membranes [T] . 
The Na+ channel transient current during a depolarizing voltage clamp may be 
described by the expression rn?h where the activation variable m and inactiva¬ 
tion variable h satisfy first order rate equations with rate functions dependent on 
the potential difference across the membrane. Support for the assumption that 
activation and inactivation are separate processes was provided by the removal 
of Na-|- inactivation from the squid axon membrane by the internal perfusion 
of pronase without affecting activation kinetics [5]. However, there is a delay 
in the onset of Nad- channel inactivation that is dependent on the time-course 
of channel activation, and Nad- channel inactivation partially immobilizes the 
gating charge associated with activation, and it was assumed that the volt¬ 
age dependence of inactivation was derived from the Nad- activation process 
[3]. There is also a delay in the recovery from inactivation that is dependent 
on the time-course of deactivation, and the rate of recovery from inactivation 
saturates for large hyperpolarizing potentials [1], and therefore, activation and 
inactivation are interdependent or coupled processes. 

The Nad- channel protein is comprised of four domains DI to DIV, each 
containing six alpha-helical segments SI to S6, and in each domain the voltage 
sensor, the S4 segment, has positively charged residues located at every third 
position. The re-entrant loop between S5 an S6 forms the ion-selective filter 
at the extracellular end of the pore, whereas the intracellular end of the pore 
is formed by the S6 segments. The inactivation gate is an IMF motif that is 
positioned on an intracellular loop between Dill and DIV, and interacts with 
and blocks the flow of ions through the inner mouth of the pore [5]. Based on 
voltage clamp fluorometry, in response to membrane depolarization, the trans¬ 
verse motion of the charged S4 segments of the Nad- channel domains DI to 
DIB is associated with activation, whereas the slower movement of DIV S4 is 
correlated with inactivation HH]. This may occur for small depolarizations 
when the ion channel is usually closed (closed-state inactivation) or for larger 
depolarizations when the S4 segments of the domains DI to D3 are activated 
(open-state inactivation). 

In a naturally occurring paramyotonia congenita mutation of the outermost 
arginine residue of DIV S4 in the human muscle Nad- channel, the inactivation 
rate is decreased with little voltage dependence for moderate depolarizations 
and therefore, the voltage dependence of inactivation is dependent on 
charged residues in the S4 segment of the DIV domain. The voltage dependence 
of the open to inactivated transition was also demonstrated by comparison of 
gating current measurement in wild-type and ApA toxin modified cardiac Nad- 
channels m- By measurement of the OFF gating charge during repolarization 
in an inactivation modified mutant of the human heart Nad- channel, it was 
estimated that the DIV S4 sensor contributes approximately 30% to the OFF 
charge, approximately 20% may be attributed to the DIB S4 sensor when the 
inactivation gate is intact, and the rate-limiting step is the motion of the DIV 
S4 sensor and not the unbinding of the inactivation gate m 
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In order to account for the effect of double-cysteine mutants of S4 gating 
charges on the ionic current of the bacterial Na-I- channel NaChBac, it has been 
proposed that at least two transitions are required during the activation of each 
voltage sensor m- This conclusion is consistent with an earlier result that cross- 
linking a DIV S4 segment from the extracellular surface inhibits inactivation 
during membrane depolarization whereas cross-linking the same segment from 
the inside inhibits activation of the Na-I- channel, and therefore, the DIV S4 
sensor translocates across the membrane in two stages m- A Nad- channel 
model that assumes that the motion of the DIV S4 sensor includes a first stage 
that is necessary for opening of the channel, and a second stage that is required 
for inactivation, provides a good description of gating and ionic currents El- 
The measurement of gating currents for charge neutralized segments in each 
domain of the Nad- channel gives additional support to the conclusion that 
the two stage activation of the DIV S4 sensor is correlated with ion channel 
inactivation m- 

In this paper, expressions for the voltage dependence of the rate of inactiva¬ 
tion and recovery from inactivation are derived by assuming that Nad- channel 
inactivation is a two stage process, where the activation of DIV S4 is correlated 
with the binding of the inactivation motif to the ion pore. However, Nad- chan¬ 
nel inactivation is, in general, dependent on the rate functions of the S4 sensors 
of each of the domains DI to DIV, and from a solution of the master equation 
for the activation of Nad- channel conductance by a single voltage sensor that 
is coupled to a two-stage inactivation process, it is shown that the voltage de¬ 
pendence of the rate functions for inactivation and recovery from inactivation 
have a similar form to empirical expressions for Nad- channels [Dill, and in par¬ 
ticular, the exponential variation exhibits saturation for both depolarized and 
hyperpolarized clamp potentials. 

INDEPENDENT ACTIVATION AND INACTIVATION OF A 
Na CHANNEL 

Inactivation of a Nad- channel may be described as the transverse motion of 
the charged S4 segment of the domain DIV, with rate functions and /3i, fol¬ 
lowed by the binding of an intracellular hydrophobic motif which blocks conduc¬ 
tion through the ion pore, with rate functions % and 6i (see Fig. 1). Assuming 
that the transition of the DIV S4 segment across two potential barriers occurs 
within an energy landscape, it may be shown from a solution of the Smolu- 
chowski equation [mull HI] that the occupation probabilities of the permissive 
states hi, h 2 and the inactivated state /13 are determined by 

= -aihi{t) + I3^h2{t), (I) 

= a^hl{t) + Sihsit) - {Pi+ 'yi)h2{t), (2) 

= lih2{t) - 6ih3{t), (3) 

where each rate is an exponential function of the membrane potential. If the 


dhi 

dt 

dh2 

dt 

dhs 

dt 
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Na+ channel is depolarized to a clamp potential V from a large hyperpolarized 
holding potential (hi( 0 ) = 1, and /i 2 ( 0 ) = ^ 3 ( 0 ) = 0 ), and it is assumed that 
the first forward and backward transitions are rate-limiting [HI HO] (A > (5* 
and 7 i » a^), the solution of Eqs. (P) to (|31) is jH] 




aili 

UJ 1 OJ 2 


acii 


— UJ2) 


exp (—Wit) — 




UJ2iuJl — UJ2) 


exp(—u}2t). 


( 4 ) 


where Wi -I- Si{ai + Pi))/(ji + Pi) and a ;2 sa 7i -f A > • 

However, if the DIV S4 sensor is initially in the inactivated state {hi{0) = 
^ 2 ( 0 ) = 0, and /i 3 ( 0 ) = 1), and the membrane is hyperpolarized to a potential 
H, it may be shown that 


hsoit) = 


ctiji -I- 5i{ai + Pi) exp(-wit) 


Q!i7i -f 5i{ai + Pi) 
Eqs. dH) and ([5|) are solutions of the rate equation P 

dhs 


dt 


— Ph,2 — {0ih,2 + Ph,2)^3, 


where 


<y-h,2{y) 

ph,2iy) 


dpOii “t“ Pi 

1% + Pi 

OLi'^i 


h + Pi 

Therefore, the probability of the permissive states h = hi + h 2 satisfies 


( 5 ) 

( 6 ) 

( 7 ) 

( 8 ) 


777 — 2 “ ((^h,2 + Ph,2)h, (9) 

at 

and Eqs. © and provide a good fit to the empirical inactivation rate 
functions at and Ph for the squid axon Na-f channel P (see Eig. 2). 

Based on the measurement of a rising phase of the gating current in a squid 
axon membrane and the chemical structure of a Na channel, at least two tran¬ 
sitions are required for the activation of each voltage sensor [121 [22]. Therefore, 
for each voltage sensor from domains DI to Dill, assuming no cooperativity 
between sensors, the occupation probabilities of the closed states mi, m 2 and 
the open state m (see Fig. 4) are determined by 


dm I 
dt 
dm 2 
dt 
dm 
dt 


-aami{t) + Pam2{t), 

aami{t) +5am{t) - {Pa +^a)m2{t), 

lam2{t) - 5 am{t). 


( 10 ) 

( 11 ) 

( 12 ) 


where the transition rates aa, Pa-, la and 5 a are exponential functions of the 
membrane voltage V. If we assume that Pa ^ 6a and la ^ aa [Hnn], from the 
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solution of Eqs. m to m during activation ( toi ( 0 ) = 1), and deactivation 
(to( 0) = 1), it may be shown that the open state m may be approximated by 


rriAit) 

moit) 


Otala 


-[1 - exp(-a;it)], 


o-ala + 5a{aa + Pa) exp(-a;it) 
~b ^a(cfco “t“ /3a) 


(13) 

(14) 


where the low frequency wi ~ (7aQ!a + 5a{oia + Pa))/{la + Pa)- Eqs. (I13p and 
(fni) satisfy the rate equation [T] 


dm 

dt 


— ^m,2 (o^m,2 “f /3m,2)^; 


and the rate functions 


(15) 


Cira,2{V) 

Pm,2{V) 




1 + Pa/la' 
da (o^a “t“ Pa) 
la + Pa 


(16) 

(17) 


provide a good fit to the empirical functions am and Pm for the squid axon Na 
channel [1] (see Eig. 4). However, the cooperativity between the S4 sensors 
in domains D1 to Dill also contributes to the voltage dependence of the Na+ 
conductance rate functions, and more recent models have adopted exponential 
functions for both am and Pm [I]- 

COUPLED MODELS OF ACTIVATION AND INACTIVATION 
OF A Na CHANNEL 

The time-dependence of the Na-I- current in the squid axon may be expressed 
as m^h where the activation variable m(t) and inactivation variable h(t) satisfy 
the rate equations [T] 

djTfi 


dt 

= ah - {ah + Ph)h. 
dt 


(19) 


The Hodgkin-Huxley (HH) description of the Na current is equivalent to an 
8-state master equation where three independent voltage sensors may activate, 
and inactivation may occur from the open state or from each of the closed states 
[23]. In this section, we assume that the activation of a single voltage sensor 
regulating the Na channel conductance is coupled to a two-stage inactivation 
process (see Fig. 5), and therefore, the kinetics may be described by a master 
equation where the occupation probabilities of the closed states Ci and Ai , the 
open states O and A 2 , and the inactivated (or blocked) states Bi and B 2 are 
determined by 


dC^ 

—— = —(ai-l-Q!o)C'i(t)-l-/3oO(t)-|-/3iAi(t) 

dt 


(20) 






— 0 :oCl{t) — (Po + Ct2)0{t) + /32^2(0 (21) 

= Q!iCi(i) — {cKA + /3i + 71)^1 (l) SiBi(t) + /3a^2(1) (22) 

= a20(i) - (;5 a+ /I 2 +72)^2(i)+^ 2 ^ 2 ( 1 )+aA^i(l) (23) 
at 

= 7iAi(t) - (as+(5i)Si(i)+/3BB2(i) (24) 

= 72 ^ 2(0 + asBilt) - (/3b + 62 )B 2 (t), (25) 


and the transition rates are functions of the membrane voltage V. 

Assuming that the first forward and backward transitions for inactivation 
are rate limiting, j3j 2> Sj and 2> aj, for j = 1, 2, Ai and A 2 satisfy 


Ai 

A 2 


ttiCi + 


,Si + 7i 
a20 + S 2 B 2 
/32 + 72 


(26) 

(27) 


and therefore, Eqs. (I20|) to (j25p iii3,y bo reduced to & four sts-te muster eQuution 
(see Fig. 6) 


—— — —{pi + (^o)Ci{t) (28) 

at 

^ = aoCi{t) - (/3o + P2)0{t) + cr2^2(t) (29) 

dB 

= piC'i(t) - (as + cri)Bi(t) + PBB 2 {t) (30) 

dB 

= P20{t) + asBiit) - {/3 b + (J2)B2{t), (31) 


where the derived inactivation rate functions 


Pi 




ai7i 0272 

/3i + 7i ’ /32 + 72 ’ 

(5i(/3i+Q;i) (32(/32+Q!2) 

-,(J2 ~ - 

,di + 7i /32 + 72 


(32) 

(33) 


If Na+ activation and inactivation are independent processes (oq = ^b, 
/ 3 o = I^Bj Pi = P2 = P and CTi = CT 2 = cr), and the Na+ channel is depolarized 
to a clamp potential V from a large hyperpolarized holding potential, the open 
state 0{t) = m{t)h{t) where the activation and inactivation variables 


m{t) = -—— (l-exp[-(ao + /3o)i]), (34) 

oo + fJO 

hit) = ■’ + i>^M-iP+At) ^ ( 35 , 
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and satisfy the rate equations (fT51) and (1191) when ao = ctm, / 3 o = Pm, P = Ph 
and a = ah (see Fig. 7). Similarly, if the Na+ channel is hyperpolarized 
to a clamp potential V from a large depolarized holding potential, Ci(t) = 
(1 — m(t))h{t) (see Fig. 8 ) where 


m{t) 

h{t) 


o-o + Po exp(-(ao + Po)t) 
ao + Po 

—^(1 -exp[-(p + cr)t]). 

p + a 


(36) 

(37) 


In the HH model, the inactivation rate functions are assumed to be indepen¬ 
dent of the Nad- conductance activation functions but, by contrast, the rate of 
recovery for inactivation in hippocampal neurons is also dependent on the Na-f 
channel deactivation functions [4]. 

More generally, the Nad- current during deactivation of the channel is very 
small [3], and therefore, the deinactivation rate a 2 ~ 0. For a depolarizing 
clamp potential V from a large hyperpolarized holding potential, the solution 
of Eqs. (051) to (OTl) when pi = p 2 = p is 


Ci{t) = Cis+ - Po - P2)exp{-u}jt) (38) 

0{t) = Os—'^^^ikj+iaoexp{—ujjt) (39) 

Bi{t) = Bis+ T,^^-i^kj+iXi{ujj)exp{-u}jt) (40) 

B2{t) = B2s + ^^=ikj+iX2{ujj)exp{-ujjt), (41) 


where Cu = kiaiPsiPo + p), Og = kiaiPsao, Bis = kipPsiao + Po + p), 
B 2 s = kipasiao + Po + p) + paocri, = [i.ao + Po + p){(^iPb + p{aB + 
Pb)) + pcriao], 


k2 

h 


ki 


Xiico) 

X2(uj} 


1 - kiaiPB0J2 d- k 4 {uJ 2 - UJ3) 

UJi - UJ2 

1 - kiaiPBiOi d- ^4(1^1 - 
UJ2 — Wl 

_ -Bis{ui2 - c^i) d- r2Xi{uJi) - riXi{uj2) _ 

(a ;2 — a;i)Xi(a; 3 ) d- (ws — ijJ2)Xi{u}i) — (wa — u}i)Xi(u}2) ’ 
-paoPB - pjuJ - Pb){'^ - Pq- p) 

- u:{aB + Pb + cri) + aiPB 
-paBjap + Po + p) - papcTi d- puj{ao + ob) 

- uj{aB + Pb + cri) + ctiPb 


(42) 

(43) 

(44) 

(45) 

(46) 


For depolarization clamp potentials, a ;2 ~ ao + Po + P, ^3 ~ aB + Pb + <^ 1 , and 


~ ah -\- Ph, 


(47) 


where the inactivation rate 


Ph = 


_ P _ 

1 d- (Jil{aB d- Pb) 


p<J\ao 

{ao + Po + p){aB + Pb + cti ) ’ 


(48) 
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(49) 


and the rate of recovery from inactivation 

PbUi 

“'i =-Tfl—^-• 

01B + PB + O'! 

The second term in Ea. (l48p only makes a contribution to Ph for small clamp 
potentials. If ai = 02 and are independent of V, uji saturates for both large 
positive and negative clamp potentials (see Fig. 9). From Fq. (l3^ . as ^4 ~ 0 
for a depolarizing potential, we may write 


0{t) 


ap / 0^(1 - exp[-(Qo + Po+ p)t]) ^ 

<^0 + Po V Oih + Ph 

Ph exp[-(Q!ft + Ph)t]{l - exp[-(ao + Po)t]) 

0!h + Ph 


(50) 

(51) 


and therefore, 0{t) 
and 


m{t)h{t) (see Fig. 10) where m{t) is defined in Fq. (IM)) 


h{t) 


ah + Ph exp(-(Q^ + Ph)t) 
ah + Ph 


(52) 


That is, the HH description of the Na+ current in terms of the variables m{t) 
and h{t) is an approximation to an expression that may be derived from a 
solution to a coupled model of Na+ activation and two stage inactivation for 
which the the deinactivation rate tT 2 ~ 0 . 

For a moderate hyperpolarizing clamp potential from a depolarized holding 
potential, the inactivation rates pi, P 2 ~ 0, and the solution of Fqs. (1^ to m 
is (see Fig. 11) 


Ci{t) 


Oit) 


i?2(t) 


Po 


exp(-..it)-ii^exp(-a; 2 t) 


ao + Po uji — UJ 2 ^ ' ' ijJi — UI 2 

aotiJi(jJ2 + i^30'2{Po — aB — cri) 


UJ 3 {uJl - UJ3){UJ2 - UJ 3 ) 

ap Y 2 {uji) 


■ exp(-w 3 t) 


exp(—Wit)- exp{—uj 2 t) — 


ao + Po U!i — UJ 2 ‘ ' ' wi — a ;2 

aoWiW2 + ljJ3^2{Po — aB — CTi) 


a; 3 (a;i - uj 3 ){ijJ 2 - W 3 ) 


Pb 


UJi — UJ 2 
wi — Ob — ui 
Wl — W 2 


■ exp(—Wit) + 


Pb 


■exp(-w 3 t), 
■ exp(-w 2 t). 


Wi — W 2 

exp(—Wit) — ———— exp(—W 2 t), 

Wl — W 2 


(53) 


(54) 

(55) 

(56) 


where W3 = ao + Po , and wi, W2 are solutions of the characteristic equation 

w^ - w(aB + /3 b + CTi + 0 - 2 ) + ctiPb + cr2(aB + cti) = 0, (57) 


-P0UJ1UJ2 + (^{ctiPb + CTlPo) 
w(w - ao - Po) 

WiW 2 (ct 2 — ao) — WtT2(/3B + 0-2 — tto) 
w(w - ao - Po) 


(58) 

(59) 
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If Pi = P 2 , from the application of microscopic reversibility to the four state 
system (see Fig. 6 ), /3b = / 3 oQ!B 0 ' 2 /ao 0 'i ^ Po for cti (T 2 ~ 0. From 
Eq. (I57)) . assuming that as “C /3 b ^ ui for a small hyperpolarizing potential, 
the lowest frequency uji ~ /3b, whereas for /3b ^ cti, ~ cti (see Fig. 9). 
The conclusion that, for small hyperpolarizing potentials, the recovery rate for 
inactivation a/j ~ /3 b oc /3m is supported by the HH data where ahiV) and 
(3m{V) have a similar voltage dependence and f3m{V) ~ 57ahiV) [1]. However, 
if as is an exponential function of V such that ob + /3b 3> cti [1] (see Fig. 12), 
from Eq. (153, W 2 ~ ttB + /3b + CTi, and 


Wl 


CTiPb 

aB + /3 b + CTi ‘ 


(60) 


Therefore, the voltage dependence of the rate of recovery from inactivation 
is determined by the deinactivation rate Ui and the Na+ conductance de¬ 
activation functions [J]. For small hyperpolarizing potentials (/3b 'C ctb), 
oJi ~ (Jij3B/[oLB + ai), and may be approximated by an exponential function of 
V [1] but saturates at a more negative potential when /3 b ^ cti ub (see Fig. 
12 ). 

From Eq. ()53p . we may write 


Ci{t) 


Po 

Oio + Po 


^1 — exp(—Wit) 


^ wi(l - exp[-(a ;2 

W2 — Wi 



(61) 


and therefore, dCi/dt{0) = 0 and there is a delay in the recovery from inacti¬ 
vation mm- However, for large negative potentials, W 2 ~ /3b ^ wi « cti, and 
Eq. (|6T|l reduces to the HH expression 


Ci{t) ^ [1 - ms]h{t), (62) 

where mg = aoj/ao /3o) and /i(t) = 1 — exp(— wit). 

CONCLUSION 

Hodgkin and Huxley described the voltage dependence of the Na-|- channel 
inactivation rate and the rate of recovery by exponential functions which for 
the inactivation rate saturates for a moderate depolarizing potential [I]. Based 
on the absence of a gating current that corresponds to the time course of the 
inactivated Na-I- current, it was assumed that the transition rates governing 
inactivation were voltage-independent, and that the macroscopic voltage de¬ 
pendence of inactivation derived from the Na-I- channel activation process [3]. 
However, the voltage dependence of the inactivation rate is dependent on the 
charge on the S4 segment residues of the DIV domain mini [in], and gener¬ 
ally, only has a minor contribution from the Na-I- conductance activation rate 
functions. The voltage dependence of the rate of recovery from inactivation 
saturates for a large hyperpolarizing potential, and has been attributed to the 
Na-I- channel deactivation rate functions |4]. 

In this paper, assuming that Na-I- channel inactivation is a two stage process, 
where the activation of DIV S4 is correlated with the binding of the inactivation 
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motif to the ion pore, we show that during a voltage clamp of the Na+ channel, 
the solution of the master equation for the inactivation process may be approx¬ 
imated by the solution of a rate equation. The backward transition rate is, in 
general, an exponential function of the membrane potential V, and the forward 
rate may be expressed as an exponential function for small depolarizations but 
approaches a saturated value for a larger clamp potential, reflecting the voltage 
independence of the rate limiting step. 

If the processes of Na activation and inactivation are independent, and the 
activation of a single voltage sensor that regulates the Na channel conductance 
is coupled to a two-stage inactivation process, the open state probability 0(t) 
during a depolarizing clamp potential, derived from an analytical solution of 
a four state master equation, may be expressed as m{t)h(t) where m(t) and 
h{t) satisfy rate equations for activation and inactivation, and the voltage de¬ 
pendence of the rate of inactivation provides a good approximation to the HH 
function Ph [I]. However, the voltage dependence of the rate of recovery from 
inactivation is dependent on the rate functions for the DIV sensor and not the 
Na-|- channel conductance deactivation functions, as observed experimentally in 
hippocampal neurons [4]. 

Based on the measurement of a very small ion channel current during de¬ 
activation [?], the Na-|- channel deactivates before recovery from inactivation, 
and hence the deinactivation rate (J 2 to the open state is smaller than the cor¬ 
responding rate cri from a deactivated state . A more general expression for the 
voltage dependence of the rate of recovery from inactivation may be determined 
that is approximated by an exponential function for physiological potentials [1] , 
but for more negative hyperpolarizing potentials, it approaches a limiting value 
equal to the deinactivation rate cti. The inactivation rate is, in general, also 
dependent on the Na-|- conductance activation rate functions, as well as the rate 
functions for the DIV S4 sensor, but for a moderate depolarizing potential it 
reduces to a two-stage expression when the DIV S4 activation rate pi = P 2 - 

We conclude that general expressions for the voltage dependence of rate 
functions for inactivation and recovery from inactivation may be determined 
from a master equation for interdependent Na-|- channel activation and inac¬ 
tivation, and are consistent with the empirical data from the squid axon [T], 
hippocampal neurons [4] and Navi.4 ion channels [15]. When Na-I- conductance 
activation is regulated by a single voltage sensor, the HH description of the Na-|- 
current in terms of the activation and inactivation variables m{t) and h(t), is an 
approximation to an expression derived from the solution to a coupled model of 
Na-|- activation and two-stage inactivation where the deinactivation rate to the 
open state (72 ~ 0. 
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«! Yi 

hi ^ h2 ^ ha 
/3i (5i 

Figure 1: Inactivation model of a Na+ ion channel, where the occupation 
probabilities of the permissive states hi{t), /i 2 (t), and the inactivated state 
hsit) satisfy a master equation, and the rate functions ai, (3i, ji and Si, are 
voltage-dependent rate functions between states. 
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Figure 2: The derived rate functions ah 2 (V) and (3h 2 (V) (solid line) in Eqs. 
© and ([5]) provide a good approximation to the HH rate functions (ms 
ah = 0.07exp(—(E-|-50)/20) and j3h = 1/(1 -|-exp(—(20-1-h)/10)) (dotted line) 
when the inactivation rate functions are ai(V) = 1 <C 7 i(h) = exp(3), and 
l3i{V) = exp(-2.5(l/ - 10)/25) > d^(y) = 0.07exp(-(E -h 50)/20). 
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Figure 3: Activation model of a Na+ channel, where the occupation probabili¬ 
ties of the closed states mi, m 2 , and the open state m satisfy a master equation, 
and Q!q, Pa, 7a and Sa are voltage-dependent rate functions between states. 
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Figure 4: The derived rate functions am 2 {V) and /3m 2 iV) (solid line) in Eqs. 
CH) and (El) provide a good approximation to the HH rate functions (ms 
Qfm = 0.1(l/+25)/(l+exp[—0.1(14+25)]) and /3m = 4exp[—(E + 50)/18] (dotted 
line) when the activation rate functions are aa{V) = 2.7exp(0.27(14+50)/25) <C 
j^(V) = 30exp(0.25(V + 50)/25), and Pa{V) = 251 exp(-0.95(14 + 50)/25) > 
5„(14) = 4.65exp(-0.9(E + 50)/18). 
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Figure 5: Six state system that describes Na+ conductance activation between 
states Cl and O, Ai and A 2 , and Bi and B 2 , is coupled to a two-stage Na-|- 
inactivation process between states Ci and Bi, and O and B 2 . 
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Figure 6 : The six state system of Fig. 5 may be approximated by a four state 
system when /3j ^ and 7 ^ ^ , for j = 1,2 where pj and Uj are derived 

rate functions for a two-stage Na-I- inactivation process, defined in Eqs. (IMI) 
and (l33)) . 
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Figure 7: For a four state system where activation and inactivation are inde¬ 
pendent, the open state probability 0{t) (solid line) during a depolarizing clamp 
potential, is equal to m{t)h{t) (dashed, dotted or dot-dashed line) where m{t) 
and h(t) are solutions of rate equations for activation and inactivation, and ao = 
= 0.1{V + 25)/(l - exp[-CF + 25)/10]), Po = Pm = 4exp[-(y + 50)/18], 
p = = 1/(1 -p exp[-(20 -b ld)/10]) and a = an = 0.07exp[-(y -b 50)/20] 

(ms-i) p. 



Figure 8: For a four state system where activation and inactivation are in¬ 
dependent, the closed state probability Ci{t) (solid line) during a hyperpo- 
larizing clamp potential, is equal to [1 — m{t)]h{t) (dashed, dotted or dot- 
dashed line) where 7 n{t) and h(t) are solutions of rate equations for activa¬ 
tion and inactivation, and ao = ctm = 0.1(14 -b 25)/(l — exp[— (14 -b 25)/10]), 
Po = Pm = 4exp[-(14 -b 50)/18], p = Ph = 1/{1 + exp(-(20 -b 14)/10)) and 
a = ah = 0.07 exp(—(F -b 50)/20) (ms“^). 
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Figure 9: Voltage dependence of the HH Na+ channel inactivation rate func¬ 
tions Ph = 1/(1 -I- exp(—(20 -I- V)/10)) and ah = 0.07exp(—(V -|- 50)/20) 
(dashed line) may be approximated by analytical expressions in Eqs. (I48p 
and (HHl) (solid line) derived from a master equation for a four state sys¬ 
tem where activation and two stage inactivation are interdependent, and by 
the voltage dependence of the lowest frequency of the system determined nu¬ 
merically (dotted line) where the rate functions are ai{V) = a 2 {V) = 1, 
MV) = MV) = exp(3), P,{V) = P 2 {V) = exp(-2(V - 10)/25) <5i(V) = 4.2, 
52{V) = 0, ao = 0.1(V -b 25)/(l - exp[-(V -b 25)/10]) = as/S and po = 
4exp[—(V -b50)/18] = 83/3 b (ms“^). 
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Figure 10: For a four state system where activation and inactivation are 

interdependent, during a depolarizing clamp potential, the open state proba¬ 
bility 0{t) (solid line) may be approximated by m{t)h{t) (dashed, dotted or 
dot-dashed line) where m{t) and h{t) are solutions of rate equations for acti¬ 
vation and inactivation, and ai{V) = a 2 {y) = 1, 71 (F) = 72 (F) = exp(3), 
/3i(tv) = / 32 (F) = exp(-2(F - 10)/25) di(F) = 4.2, 52{V) = 0, ao = 0.1(F -h 
25)/(l - exp[-(F -h 25)/10]) = as/i and /3o = 4exp[-(F -h 50)/18] = 83/3ij 
(ms“^). 
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Figure 11: For a hyperpolarizing clamp potential of a four state system where 
activation and inactivation are interdependent, the closed state variable C'i(t) 
(solid line) may be approximated by [1 — ms\h{t) (dashed, dot or dot-dashed) 
where rus = aol{oco + Po) and h{t) is a solution of a rate equation for inac¬ 
tivation, and rate functions are ai(F) = a 2 {V) = 1, 71 (F) = 72 (F) = exp(3), 
/3i(F) = / 32 (F) = exp(-2(F - 10)/25) di(F) = 4.2, (52(F) = 0, ao = 0.1(F -h 
25)/(l - exp[-(F -h 25)/10]) = as/S and /3o = 4exp[-(F -h 50)/18] = 83/3ij 
(ms“^). 
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Figure 12: Voltage dependence of the HH Na+ channel inactivation rate func¬ 
tions Ph = l/(l-l-exp(—(20-|-V)/10)) and an = 0.07exp(—(V-|-50)/20) (dashed 
line) may be approximated by analytical expressions (solid line) derived from 
a master equation for a four state system where activation and two stage in¬ 
activation are interdependent, and by the voltage dependence of the lowest 
frequency of the system determined numerically (dotted line) where the rate 
functions are ai(V) = a 2 {V) = 1, 7 i(V) = 72 (V) = exp(3), Pi{V) = P 2 {V) = 
exp(-2.5(V - 10)/25) di(V) = 4.2, 52{V) = 0, ao = as = 5.4exp(-0.3V/25) 
and Po = 0.5exp[—V/18] = 98Pb (ms“^). 
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